Chaos and Interacting Electrons in Ballistic Quantum Dots 
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We show that the classical dynamics of independent particles can determine the quantum properties 
of interacting electrons in the ballistic regime. This connection is established using diagrammatic 
perturbation theory and semiclassical finite-temperature Green functions. Specifically, the orbital 
magnetism is greatly enhanced over the Landau susceptibility by the combined effects of interac- 
tions and finite size. The presence of families of periodic orbits in regular systems makes their 
susceptibility parametrically larger than that of chaotic systems, a difference which emerges from 
correlation terms. 



The connection between classical dynamics and wave 
interference has attracted attention in many fields of 
physics recently including atomic, mesoscopic, and 
optical physics. A central question is to what extent 
the quantum properties of classically regular and chaotic 
systems differ. On the whole, this question has been ad- 
dressed for non-interacting systems. It is now known that 
many quantum properties are in fact strongly influenced 
by the nature of the classical dynamics- the density of 
states, the quantum corrections to the conductance, and 
the optical absorption to name a few. 

We wish to address this question for interacting sys- 
tems and, in particular, to investigate the role of the clas- 
sical dynamics of the non-interacting system in this con- 
text. If the interactions are strong, the non-interacting 
classical dynamics will be of little relevance. However, 
if the interactions are short-range and not too strong, 
the non-interacting classical dynamics may be impor- 
tant, and its role can be assessed with perturbation the- 
ory. This regime is physically relevant: it applies to a 
high density two-dimensional electron gas in which the 
quasi-particles interact weakly through the short-range 
screened Coulomb interaction. We find that at first or- 
der in the interaction there is a difference between reg- 
ular and chaotic systems, but one which is only numer- 
ical, not qualitative. Intriguingly, as the perturbation 
theory is carried out to higher orders a qualitative differ- 
ence emerges: thermodynamic properties scale differently 
with Fermi energy for chaotic and regular systems. This 
correlation effect shows that the nature of the classical 
dynamics can have a substantial effect on the quantum 
properties of an interacting system. 

To be specific, we study the magnetic response of an 
ensemble of ballistic quantum dots formed from a two- 
dimensional electron gas. Recent fabrication progress has 
made possible phase-coherent electronic microstructures 
much smaller than the mean free path. In these "bal- 



listic" quantum dots, one can think of electrons mov- 
ing along straight lines between specular reflections off 
the confining potential. Because this motion is qualita- 
tively different from that taking place in bulk materials, 
a variety of new behavior has been observed experimen- 
tally H]. In particular, the magnetic susceptibility of 
an ensemble of ballistic squares has been measured 
and a large enhancement over the Landau response X-t 
was found. First attempts to understand this experiment 
within non-interacting models pointed to the importance 
of the classical dynamics [||J4[|^ . The inclusion of interac- 
tions in such systems is our main concern, though much 
of the discussion applies to ballistic structures in general. 

For the magnetic response, the high-density expansion 
(RPA) of the thermodynamic potential ||^ has to be ex- 
tended by including Cooper-like correlations, as carried 
out previously for disordered metals |0,||,D • Such expan- 
sions are typically used beyond the high-density limit and 
yield reliable results for the bulk provided some sets of 
terms are properly resummed. We continue to follow 
this point of view for quantum dots, where the "small 
parameter" Ts = ro/ao is about 2. (Trrg is the average 
area per electron, and ag is the Bohr radius in the ma- 
terial.) We show that these expansions are particularly 
insightful when combined with a semiclassical approxi- 
mation from which the connection to the nature of the 
classical dynamics can be made. Thus, we will assume 
that kpa^ 1 (a is the size of the microstructures and kp 
the Fermi wavevector) and that the magnetic field B is 
classically weak (cyclotron radius ^ a). 

Semiclassical approach.- The perturbation expansion 
for the interaction contribution to the thermo- 
dynamic potential Q yields the magnetic susceptibility 
through X = {—l/a'^)d^^l/dB^. The dominant terms 
are shown in Fig. 1. The screened Coulomb interac- 
tion (wavy lines) is treated as local |10|, U{r — r') = 
\oN{0y^5{r - r'), with N{0) the density of states and 
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FIG. 1. Leading Cooper-channel diagrams for the interac- 
tion contribution to the thermodynamic potential. 

Aq = 1 identifying the order of perturbation. Straight 
lines represent the "free" finite-temperature Green func- 
tion in the presence of the confining potential, 

gr,r'ie„)^eie^)G^:^,{E^ + ie^) + e{-e,,)Gl^,{E^+ie,,) . 

Here, Ep is the Fermi energy, e„ = (2n -I- 1)tt/(3 the 
Matsubara frequencies, and G^'^ the retarded, advanced 
Green functions related by G^y{E) = [G^, ^^{E*)]* . 
Semiclassically, G^ is the sum of the contributions 
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of each classical trajectory j from r to r' [Q: 



G?AE) 



. y D.e''^^/^-''"''''^ (1) 



where Sj— p • dr is the classical action of trajectory 
{xx')~^ \d'^Sj/dydy'\ the classical density, and 



J, 



Uj a Maslov index. Using {dSj / dE)=tj and {dSj/dB)= 
{e/c)Aj, where tj and Aj are the traversal time and area, 
one finds that the contribution of trajectory j to G^ is 



G^;AE,+ten,B) = G';;^,{E,,B = Q) 

X exp [—entj/K\ x exp [i27rBAj/(j)o] (2) 

where 0o = hc/e is the flux quantum. Note that tempera- 
ture introduces time and length scales tx = Lt /vp = HP/tt 
which exponentially suppress the contributions of long 
paths through the term €ntj /h—{2n+l)tj/tT- {vp is the 
Fermi velocity of a billiard.) This provides a complete 
description in the semiclassical perturbative regime. 

We start with the first-order (Hartree-Fock) term in 
the diagrammatic expansion 



(3) 



where the trace refers to the spatial arguments of the 
particle-particle propagator M 
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{lu ^ LOn — ifi-K/P). The short-length (high-frequency) 
behavior is incorporated in the screened interaction, thus 
requiring a cutoff of the frequency sums at Ep |^ . Semi- 
classically, S]r,r' is a sum over pairs of trajectories join- 
ing r to r'. However, most pairs yield highly oscillat- 
ing contributions which, after the spatial integrations. 





FIG. 2. Typical pairs of real-space trajectories that con- 
tribute to the average susceptibility to first order in the in- 
teraction in the diagonal channel (left) and the non-diagonal 
channel (right). 

give higher order terms in l/kpU. To leading order, only 
those pairs contribute to the susceptibility whose dynam- 
ical phases ex-p[iSj{B = 0)/h] cancel while retaining a 
magnetic-field dependence. This is achieved by pairing 
each orbit j with its time reverse. The trace in Eq. (|^) 
yields a sum over closed but not necessarily periodic tra- 
jectories (see Fig. 2, left, for a square). This "diagonal" 
or "Cooper channel" is present independent of the nature 
of the classical dynamics , and we will return to it be- 
low. We first turn to an additional contribution present 
for integrable systems which is central to this paper. 

Non-diagonal channel - In integrable systems, periodic 
orbits come in families within which the action integral 
is constant. If, as is generally the case, in some part of 
configuration space two orbits of the same family cross 
at a given point, it is possible to cancel the dynamical 
phases by pairing them (Fig. 2, right). This non-diagonal 
first-order contribution involves a term for each family 
of periodic orbits. For the square billiard, and at not 
too low temperature {Lt ^ 2a), only the shortest of 
these periodic orbits has to be taken into account, namely 
the family (11) with length Lu — 2\/2a shown in Fig. 2 
(right). In this case, we find for the susceptibility 
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the temperature dependence is governed by the function 
R{x) = a;/sinh(x) and the field dependence by C{ip) — 
(2(p)^-^/^ [cos(7r(/3)C(y^7r^) -I- sin(7r(^)S(y^7n^)] , with — 
Ba'^/<j)Q, and C and S Fresnel functions. As in the non- 
interacting case [^,|| , the contribution Eq. (||) is linear in 
kpa and has a temperature scale related to the length of 
the periodic orbit. Quantitatively, the non-diagonal con- 
tribution of the family (11) and its repetitions is shown as 
the dashed curve in Fig. 3. Thus the existence of a family 
of periodic orbits- a characteristic of the non-interacting 
classical dynamics- is associated with an additional first- 
order interaction contribution to the susceptibility. 

Higher-order terms in perturbation theory also contain 
non-diagonal contributions. However, in these terms the 
location of the additional interaction points is severely 
limited: they must lie on both periodic orbits to cancel 
the dynamical phases and so must be near the intersec- 
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FIG. 3. Temperature dependence of the zero-field suscepti- 
bility (solid line) for an ensemble of squares at fe_Fa=50. The 
contribution of the non-diagonal channel (dashed, family (11) 
and repetitions) exceeds that of the diagonal Cooper channel 
(dotted) at low temperatures {ksTo = Tivp /2-Ka). 

tions of the two orbits. Further analysis shows that these 
contributions are therefore smaller by a factor of l/kpa. 
By contrast, we will now show that the diagonal contri- 
bution is strongly renormalized by higher-order terms. 

Diagonal Cooper channel.- The first-order diagonal 
channel gives a contribution to x with the same depen- 
dence on kpa and T as in Eq. (^. So, to first order in the 
perturbation the difference between chaotic systems- for 
which there is only the diagonal contribution- and regu- 
lar ones- for which the non-diagonal contribution is also 
present- is numerical but not qualitative. 

However, higher-order diagrams are essential in the di- 
agonal Cooper channel, as known from the theory of su- 
perconductivity 1^,0. One should sum all terms which 
(i) do not vanish upon ensemble averaging, (ii) depend 
on B, and (iii) are of leading order in ?i ^ l/kpa. This 
yields the Cooper series shown in Fig. 1. For instance, 
(iii) is checked by h power counting, since a pair of Green 
functions scales as N{0)/'h, interactions as [iV(0)]^^, and 
Matsubara sums as h. Indeed, all terms in the series are 
of order h despite the formal expansion in Aq. Summing 
the series yields for the diagonal contribution M 



(6) 



The diagonal part S^^-* of S is a sum over over all trajec- 
tories longer than the cutoff Aq = Xp/tt [associated with 
the upper bound Ep on the Matsubara sum in Eq. (0)] : 



Lj>Ao 



R{2t,ltT) 



X exp [iATTBAj/(j)^)\ x exp[— wtj7?i] . (7) 

While we cannot diagonalize S'^^ analytically, it has the 
nice property that (except for Aq) all variations occur 



on classical scales: rapid oscillations on the scale of \p 
have been washed out. The original quantum problem 
becomes much simpler, involving only the "classical" op- 
erator S^^^; hence, we can discretize E^^' with mesh size 
larger than \p and then compute f2*^^) numerically. 

We have performed this computation for the square 
billiard, obtaining the dotted curve of Fig. 3 for x(T). 
In this curve, we can distinguish three regimes. At low- 
temperature x*-^-* is paramagnetic and decays on a scale 
similar to the non-diagonal contribution (dashed curve) , 
but has a significantly smaller amplitude. In the inter- 
mediate range, x''^'' is small and diamagnetic. Finally, 
at high temperatures x'^^ is again paramagnetic, but 
very small. This is naturally understood by associat- 
ing each regime with an order in the perturbation series. 
The low-temperature part corresponds to the first-order 
term (orbits of the type in Fig. 2, left) which is expo- 
nentially suppressed by the temperature factor R when 
Lt becomes smaller than the shortest closed orbit. At 
this point the second-order term, due to closed paths 
of two trajectories connected by interactions, takes over. 
There is no minimum length of these paths, and hence 
the second-order term is less rapidly suppressed by tem- 
perature. For repulsive interactions, the sign is opposite 
of the first-order term, thus the sign change in x^^-* • At 
even higher temperatures once Lt <IC a, this term is a 
surface contribution and the third-order term takes over. 
The latter is a bulk contribution [Q| since with three in- 
teractions flux can be enclosed without bouncing off the 
boundary. 

Renormalization scheme. - This interpretation of Fig. 3 
must, however, be reconsidered. First, as noted above the 
first-order diagonal and non-diagonal contributions are of 
similar magnitude while the low-temperature magnitudes 
in Fig. 3 are different. Second, one observes numerically 
that the terms in the perturbation series increase in mag- 
nitude with order: one is not in the radius of convergence 
of perturbation theory but in its analytical continuation. 
Both facts contradict the above picture. 

However, the interpretation is valid once the interac- 
tion entering the diagonal contribution is replaced by a 
renormalized interaction. To see this, we propose a sim- 
ple renormalization scheme where integration over short 
trajectories yields a decreased effective coupling constant. 
To that end consider a new cutoff A larger than Aq but 
much smaller than any other characteristic length (a, Lt, 
or y^(j)o/B). For each path j joining r to r' with Lj > A, 



let E^ J., denote its contribution to '^[.^J and define 



+ \l [ dri dr2 Ei:,riSri,r.Er..r' + • ■ • • (8) 



where the integration is over Aq < [r^-i — r^l < A (with 
ro = r'). Eri,r' is defined by Eq. (^ but with the sum re- 
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stricted to "short" trajectories with lengths in the range 
[Ao,A]; is obtained from T,^ ^, by continuously de- 

forming trajectory j. To avoid the awkward In in Eq. (||), 
we introduce r= (l//3)^„Tr[l + AoS|.^K^)]"\ from 
which ri^^) can be derived through 



^0 



(9) 



Replacing S by S in F amounts to a reordering of the 
perturbation expansion of F in which short paths are 
gathered into lower-order terms. Moreover, if Lj ^ A, 
small variations in the spatial arguments do not modify 
noticeably the characteristics of T,^ . Approximating 

l/47r|ri-r'|2 valid 



by J:ly in Eq. (||) and using En, 
for short paths, we obtain 
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where the running coupling constant is defined by A (A) = 
Ao/[l + (Ao/2) ln(A/Ao)]. Therefore, these successive 
steps amount to a change of both the coupling constant 
and the cutoff (since now trajectories shorter than A must 
be excluded) without changing F; that is. 



F(Ao,Ao)-F(A,A(A)) 



(11) 



Through Eq. (^), this renormalization scheme can be ap- 
plied to and so to the average susceptibility. 

In this way, we have eliminated the last "quantum 
scale" Ao from the definition of E^^^: A can be made 
much larger than Xp while remaining smaller than all 
classical lengths. Furthermore, it is qualitatively reason- 
able that the perturbation series of fi^"^^ becomes conver- 
gent when A is of order a, since by this point the spread 
in length scales causing the divergence has been elimi- 
nated. We have checked that this is true numerically, 
although this is at the border of the range to which A 
can be pushed with a reliable quantitative answer. The 
interpretation of Fig. 3 is therefore correct, once the cou- 
pling constant Ao = 1 is replaced by the renormalized 
one A(a) ~ 2/ln(fcj^a). The smallness of A(a) explains, 
first, why the magnitude of the diagonal contribution is 
reduced below the off-diagonal one in Fig. 3 and, second, 
why the diamagnetic excursion and high-temperature tail 
are small. 

Consequently, at low T {Lt ^ shortest periodic or- 
bit), the diagonal contribution is parametrically smaller 
than the non-diagonal one by a factor 1/ ln(fcj?a) because 
higher-order correlation terms reduce only the diagonal 
contribution. Therefore, regular systems, for which there 
is a non-diagonal contribution, show a magnetic response 
logarithmically larger than generic chaotic systems, for 
which only the diagonal channel is open. For compar- 
ison, we note that the non-interacting contribution ob- 



tained previously is of the same order as this inter- 
action contribution for integrable systems but smaller for 
chaotic ones. 

To conclude, we have shown that a semiclassical treat- 
ment allows one to study the high-density perturbative 
expansion of the interaction contribution to the grand 
potential for ballistic quantum dots. This semiclassical 
approach is an efficient tool to compute quantitatively 
the magnetic response. Moreover, when combined with 
a renormalization scheme, it provides an intuitive picture 
of various features specific to the ballistic regime. The 
most striking one is that the susceptibilities of integrable 
and chaotic geometries scale differently with kpO because 
of the presence of families of periodic orbits in the for- 
mer. Another unusual property, caused by the different 
T dependence of different orders in the (renormalized) 
interaction, is that with increasing temperature the in- 
teraction contribution changes sign from paramagnetic 
to diamagnetic and then back to paramagnetic. 
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